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The unique time signature of the survival probability exactly at the exceptional point parameters 
is studied here for the hydrogen atom in strong static magnetic and electric fields. We show that 
indeed the survival probability S{t) = \{tp{0)\tp{t))\^ decays exactly as |1 — atpe"'"^^*'''' where Fep 
is associated with the decay rate at the exceptional point and a is a complex constant depending 
solely on the initial wave packet that populates exclusively the two almost degenerate states of the 
non-Hermitian Hamiltonian. This may open the possibility for a first experimental detection of 
exceptional points in a quantum system. 
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I. INTRODUCTION 

Exceptional points (EP) jj], i.e., branch point singu- 
larities in non-Hermitian physical systems, where two 
complex eigenvalues degenerate and the corresponding 
eigenstates coalesce, have shown to exhibit prominent ef- 
fects not observable in their absence. Most dramatic is 
the influence of EPs in quantum mechanics, where ef- 
fects appear which are not possible in the case of Hermi- 
tian Hamiltonians with potentials describing bound state 
spectra. Although the EPs are single points in an (at 
least) two-dimensional parameter space they influence a 
whole region of parameters and lead to unusual results 
as the permutation of eigenstates for a closed adiabatic 
loop in the parameter space or a special type of a geo- 
metric phase 0. Exceptional points have recently been 
detected in number of physical applications. Most ex- 
pamples are known for optical systems such as unstable 
lasers [3|, waveguides Q and optical resonators In 
quantum systems the existence of exceptional points has 
been proven theoretically, e.g., in atomic [1, 0| or molec- 
ular ^] spectra, in the scattering of particles at potential 
barriers [9| , in atom waves Mini, and in non-Hermitian 
Bose-HublDard models jT^.ll3|. The experimental verifi- 
cation of th eir p hysical nature was achieved in microwave 
cavities Despite this success an experimental ob- 

servation in a true quantum system is still lacking. Reso- 
nances at an exceptional poin t exhibit, however, a unique 
decay behavior 0, [ij, [l6l - [l9| and it is the purpose of this 
article to demonstrate that this can open the possibility 
for the first experimental detection of EPs in an atomic 
quantum system. 

The most fundamental quantum objects which contain 
exceptional points are atoms in static external magnetic 



and electric fields. As such they are accessible to both 
experimental and theoretical methods, and thus ideally 
suited for studying the infiuence of exceptional points on 
quantum systems. Indeed, the existence of branch points 
in the resonance spectra of the hydrogen atom in crossed 
electric and magnetic fields was found numerically 
Here the two field strengths play the role of two control- 
lable parameters necessary to set the system at the ex- 
ceptional points. In this article we want to demonstrate 
that the unique time signature of the survival probability 
exactly at the exceptional point parameters also appears 
in a detectable form in spectra of the hydrogen atom. 

The article is organized as follows. In Sec. HI] we re- 
view how in every quantum system exhibiting exceptional 
points a unique time behavior of the survival probability 
leads to an unambiguous fingerprint of the branch point 
singularity. To verify the existence of this signal in a true 
quantum system we show that it appears for the hydro- 
gen atom in crossed electric and magnetic fields in Sec. 
mil Conclusions are drawn in Sec. IIVI 
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II. FINGERPRINT OF EXCEPTIONAL POINTS 
IN THE TIME BEHAVIOR OF THE SURVIVAL 
PROBABILITY 

Let us first explain the motivation of our work. It 
has been shown theoretically in a two-dimensional model 
[l6j , in optical microspirals Q , in a non-Hermitian Bose- 
Hubbard model [l^, and in complex crystals [13, [3 
as well as experimentally for Rabi oscillations in a mi- 
crowave cavity ^19'| that when the spectrum of a non- 
Hermitian Hamiltonian has an exceptional point then 
for a broad range of initial conditions the survival prob- 
ability, S{t) — \{'ip{0)\^p{t))\'^ , decays exactly as |1 — 
atpe+^^^'^^p)*/'*, where E^p is the complex energy of 
the resonance state at the EP and a is a complex constant 
depending solely on the initial wave packet. The reso- 
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nance decay rate (inverse lifetime) is defined as Fep = 
— 2Im(£^Ep) > 0. This behavior is in clear contrast to the 
purely exponential decay far away from an EP and the 
special condition is that the initial wave packet should 
populate only the exceptional eigenstate IV'ep) and its 
complimentary state \x) as obtained from the Jordan 
chain formalism such that |V-'ep)(xI + Ix)(^p| = 1 (see 
a detailed explanation in Sec. 9.2 of Ref. [2l| where the 
closure relations for non-Hermitian Hamiltonian with an 
incomplete spectrum is discussed in detail). In this ar- 
ticle we will demonstrate this behavior of the survival 
probability for a quantum mechanical system. 

Let us give here a simple explanation for this unusual 
situation. The EP is associated with a situation where 
H{X)\4>j) = Ej\%pj) such that for A — ?> Aep two eigen- 
values degenerate, Ej — Eji — )■ (i.e., upon coalescence 
Ej = Ej' = Eep) and also the corresponding eigenstates 
coalesce tpj ~^ V'j' (up to a phase factor i [2l|, i.e., upon 
coalescence ipj = i^j' = V'ep) such that 

MIV'EP) = (1) 

where 

M = H - E^p (2) 

and Fep = — 2Im(i?Ep) is the decay rate of the system. 

Let us take as a basis set consisting of ■i/'EP and its 
complimentary state x for the initial wave packet. In this 
basis the 2x2 matrix representation of M is a matrix 
M for which = 0. Therefore 

U{t <- 0) = e"*^*/'"' ^ ^-lEEpt/h^^iMt/h 

n 

= e-'''^^'^\h^2'tMt/h). (3) 

Consequently for any initial state iplt — 0) which is a 
linear combination of ipEP and its complimentary state x 
then 

= e-'-EEPtA(i _ it{^(^t = 0)|M|V'(i = 0))/n) (4) 

and for real values of {ip{t = 0)\M\ijj{t = 0)) the survival 
probability is given by 

s{t)^\m = o)\u{t^o)m^o))\^ 

= (1 + [{^j{t = 0)|M|V'(i = 0))/h]h'^) e+^i'^f-EEi- 

(5) 

For complex values of {ipit = 0)\M\ip{t — 0)) an ad- 
ditional term linear in t is added. Surly, the quadratic 
dependence for short times is not surprising. What is 
important here is that there are no terms of order higher 
than t"^ so that the time dependence remains for all 
times. This provides a unique fingerprint proving unam- 
biguously the presence of an exceptional point since the 



power series expansion of Eq. ([3]) stopping after the lin- 
ear term requires the presence of an EP. The effect even 
remains in a larger vicinity around the branch point sin- 
gularity. Observations which are similar to exceptional 
points but not connected to true branch points as narrow 
avoided crossings for Wannier-Stark resonances [l^ are 
not sufficient. 



III. NON-EXPONENTIAL DECAY OF 
RESONANCES IN SPECTRA OF THE 
HYDROGEN ATOM 

In our study the resonances are calculated numerically 
exact by the diagonalization of a matrix representation 
of the Hamiltonian. Without relativistic corrections and 
finite nuclear mass effects ^] the Hamiltonian reads in 
atomic units 

H=^p^-^ + ^^L, + ^^\x^+y^) + fx, (6) 

where Lz is the z component of the angular momentum. 
The strengths of the electric and magnetic fields are la- 
beled / and 7, respectively. We exploit the fact that the 
parity with respect to the (z = 0)-plane is a constant of 
motion and include in all our calculations only resonances 
with even z-parity. 

To uncover the decaying unbound resonance states we 
use the complex rotation method [U [11] , for which 
the coordinates of the system r are replaced with the 
complex rotated ones re"'. For the application of the 
complex rotation method to hydrogen spectra see [1^. 
This procedure renders the resonance wave functions 
square integrable so that that they are automatically in- 
cluded in the spectrum as new discrete eigenstates with 
complex eigenvalues in the matrix representation. The 
real parts of the complex energies represent their en- 
ergies and the imaginary parts their widths (lifetimes). 
After introducing complex dilated semiparabolic coordi- 
nates [l?! the Schrodinger equation of the Hamiltonian 
([6|) assumes in a basis representation the form of a gen- 
eralized eigenvalue problem 

A{-fJ)^ = 2\b\^EB^ (7) 

with a complex symmetric matrix A(7, /) and a real 
symmetric matrix B. In this equation b is the com- 
plex dilation parameter and E the complex resonance 
energy. The eigenstates can be normalized such that 
(\I'i|B|^'j) = Sij, where the round braces indicate an 
inner product in which complex parts originating exclu- 
sively from the complex dilation parameter b are not con- 
jugated, which is the appropriate inner product for com- 
plex scaled wave functions (c-product, cf. Refs. f2TI.[25j). 
Note that this normalization does not hold at an excep- 
tional point where each of the two states is orthogonal to 
itseff [21|]. 
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We first demonstrate that the decay signal of the prob- 
ability density according to Eq. 

Sn.it) = exp [2Im(i;Ep)i] |1 - i(*o|M|'J'o)i|' , (8a) 

can be found in the quantum spectrum of the hydro- 
gen atom in a large region around the exceptional point, 
where in our case the matrix M of Eq. ([2]) is given by 



M = A(7,/)/(2|6|*)-SepB. 



(8b) 



The subscript m of S^i in Eq. (|8ap is given in order to em- 
phasize the fact that the survival probability is calculated 
here using the matrix representation of the Hamiltonian 
given by Eq. (I8b|) . This is to distinguish from the direct 
evaluation of the survival probability as discussed later 
[Eq. ([T^]. For this purpose we calculate resonance spec- 
tra for several small distances S to the exceptional point 
at (/ep, 7ep) in the space of the two field strengths, i.e., 
we use 



/ = /ep(1-(5), 7 = 7ep(1~<5)- 



(9) 



For small S the two eigenstates and belonging 
to the branch point singularity can be identified clearly. 
Note, however, that due to roundoff errors in the nu- 
merical calculations S never ever gets strictly the value 
and always (5^0. As we will show in this article even 
when (5 7^ 0, i.e., when we have two almost degenerate 
states and the spectrum is complete the fingerprint of the 
EPs (strictly obtained only exactly at the coalescence of 
two eigenstates) on the survival probability is still pro- 
nounced. 

It is well known for exceptional points that 'i'l and 5*2 
converge to the single independent eigenstate with the 
phase relation 



*EP = *i = i*2 



(10) 



for (5 — )■ [2l|. A complete basis in the corresponding 
two-dimensional subspace is spanned by 'I'ep and the 
associated vector 5*3, 



(11) 



(A(7,/)-2|6|4£;EpB)*a = *EP, 



where is essential for a decay signal of the form (|8aj) 
[l6]. Despite the convergence of 'i'l and to the same 
state it is possible to extract an adequate superposition 
of ^I^ep and ^a, viz. 



^(1 + l/VS)/2 *i + Y^(l - 1/V6)/ 



2*5 



(12) 



As we get closer to the EP when S the amplitudes 
of 4'i and ^'2 approach infinitely large values. However, 
the initial wave packet *o remains c-normalized. About 
the complex normalization rather than using the conven- 
tional scalar product for calculating the norm of a vector 
state read Ref. l2]] |. 

We choose the exceptional point labeled 8 in Table 
I of Ref. [23 at the field strengths / = 0.0002177 and 
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FIG. 1. (Color online) Survival probability for the resonances 
at the exceptional point labeled 8 in Ref. [13] with an offset 
of (5 = 10"^ Shown are the direct evaluation Sd{t) according 
to Eq. p3|l and the expected form Sm{t) of Eq. (I8ap including 
a linear term in t besides the exponential decay. Note that on 
the scale of the figure the survival probabilities calculated by 
both methods are not distinguishable. Since we use here the 
c-product rather than the regular scalar product the survival 
probability can assume values larger than one. This non- 
physical behavior results from the way we normalize 'I'o [see 
Eq. p2p ]. As we will show later this problem disappears when 
the initial wave packet is prepared by using a laser excitation 
of the field-free ground state. On the c-product read in Refs. 

MM- 



7 = 0.004604, and with the complex energy E^p — 
-0.022135-0.00006878i (aU values in atomic units). The 
survival probability for the superposition (|12p is plot- 
ted in Fig. [l]for an offset 5 — 10~^. At an exceptional 
point we expect a decay of the survival probability in 
the form ([8a|) . The corresponding numerical result is 
shown with the solid red line in Fig. [l] where we found 
{^o\M\^o) = (3^3-4.58i) x 10"''. Since we use here the 
c-product [2lL |25| rather than the regular scalar product 
the survival probability can assume values larger than 
one for the mathematical choice given in Eq. p^ . Ad- 
ditionally, we calculate the survival probability directly 
without any assumption about its shape close to an ex- 
ceptional point, i.e., we evaluate 



Sdit) = 



^(*o|-B[*,)(*^|B|*o) exp {-iE,t) 



(13) 



where we use 100 eigenvectors of the matrix diagonaliza- 
tion in the energy vicinity of the EP for the basis states 
^'j . The blue dashed line in Fig. [T] shows the results of 
the latter method. As can be seen clearly both methods 
agree very well, which proves that the description of the 
decay with the linear term in (j8ap is correct. 

Obviously the polynomial contribution influences the 
decay significantly, i.e., the unique time behavior of the 
resonances at an exceptional point is a relevant effect in 
matter waves and can unambiguously be found for atomic 
resonances. 

The typical time signal of an exceptional point is even 
present at larger distances 5. In Fig. [Ifa) we plot the di- 
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FIG. 2. (Color online) (a) Survival probability S^ii) for the 
same exceptional point as in Fig. [T] but for several different 
offsets (5 = fO~* . . . fO~^. Small deviations only appear for 
the largest 5 — fO~^. In this case also a small difference to 
the shape Sm(t) at an exceptional point becomes observable, 
(b) The modulus d = (*i -|-i*2|*i +1*2) vanishes at the 
exceptional point as is expected due to the phase relation 



rect evaluation (|13p for the same exceptional point as in 
Fig. [T]but for several different offsets b — 10~® . . . 10~^. 
For all of these distances the two vectors belonging to 
the branch point are well defined. Almost all calculations 
provide exactly the same results. Only for 8 — 10^^ we 
observe a slight difference to the other calculations. At 
this distance also the validity of the matrix representa- 
tion (j8ap including only the two components associated 
with the exceptional point breaks down, however, the dif- 
ferences are still small. The corresponding line S^[€) is 
included in the figure. For all other values of 8 shown 
we checked that the results of both methods agree com- 
pletely, which demonstrates that the structure of Eq. ([3]) , 
which is only fulfilled in the presence of an exceptional 
point, survives in a larger vicinity around the branch 
point. The signal keeps its unique structure. To verify 
that we obtain the correct vectors and \E'2 in all cal- 
culations we plot the modulus d — -f z^'2|^'i -f i^-i) 
in Fig.[2Ub). According to the phase relation ((TU)) d must 
vanish in the limit (5 — > 0. Exactly this behavior is found. 

So far we demonstrated that it is possible to find an 
adequate superposition of the two eigenvectors, however, 
we want to show furthermore that this signal can be ex- 
cited in a realistic case. Is it possible to occupy such a 
superposition in an experimental situation? To investi- 
gate this question we assume a hydrogen atom in external 
fields, where the electron is in the orbital 2p, m = 0, for 



FIG. 3. (Color online) (a) Occupation probabilities \Ai\^ in 
arbitrary units for resonances in an energy vicinity of the 
exceptional point. The two resonances connected with the 
exceptional point are marked with blue (filled) circles. On 
the scale of our plot the two almost degenerate states are not 
distinguishable, (b) Survival probability for the excited state 
"tp, which is prepared by using a laser to excite the field- free 
ground state [see Eq. (|16p ]. The division by the exponential 
part in the inset demonstrates the presence of the polynomial 
contribution. The behavior of the survival probability as pre- 
sented here is exactly as the analytical expression given in Eq. 



which any perturbation due to the fields we use can be 
ignored. The eigenstates at the exceptional point are ex- 
cited with a laser polarized linearly along the direction of 
the static magnetic field. We use a Gaussian pulse shape 
of the form 

£;(w) - exp (-cr(w - wo)^) , Wo = R.e(£;bp)-£^i, (14) 

where E\ is the energy at the initial state (2p, to = 0). 
The width was chosen to be cr = IOOO/wq. Then the 
occupation amplitude for a transition to eigenstate of 
the Hamiltonian is 



Ai = dujE{uj) 



E2p^ — Ei,p + 



(15) 



with the dipole operator D for the present choice of the 
light pulse. 

Fig.[3Ja) shows the occupation probability \Ai\'^ versus 
real part of the energy for the states in the vicinity of the 
branch point. One can see that the two states connected 
with the branch point (marked with filled blue circles) 
have an occupation probability almost three orders of 
magnitude larger than all other states. This does not tell 
us, however, whether or not an adequate superposition 
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of the two dominating states similar to the mathematical 
case in Eq. ([T^ can be achieved. Thus, we construct the 
normalized state 

*F = ^^^=, ^-F^^A,*, (16) 

^(^'f|B|*f) 

occupied by the laser with the states shown in Fig. [3^a) . 
The survival probability is calculated according to Eq. 

with instead of 4'o- Fig. |31Jb) shows the results. 
The small oscillations are due to the weaker excitations 
of the neighboring states. They disappear for a pulse 
denser in frequency space. The dominating signal is still 
formed by the two states associated with the exceptional 
point. The linear part in the time behavior (|8ap is weaker 
than in the mathematical case of Eq. (fT2)) . however, it 
is present and is expressed in the non-exponential decay. 
After the division by the exponential part the polynomial 
contribution of the physical (observable) survival proba- 
bility calculated for the initial wave packet [see Eq. 
([TO)) ] is given by 

Sd^t) - S'd(i)/exp(2Im(i;Ep)i) . (17) 

To demonstrate that the origin of this signal is in fact 
the structure (|8aP originating from an exceptional point 
we calculated the matrix element (^fIA^I^'f) = (0.226-1- 
5.25i) X 10"^ The line |1 - i(4'o|M|^'o)t|^ is not distin- 
guishable from the full numerical result presented in Fig. 

Hb). 



IV. CONCLUSION 

We proved in this article that any quantum system ex- 
hibiting exceptional points shows a time evolution of the 
form ^ for two resonances exactly at the EP, i.e., the 
decay includes a quadratic term as in Eq. ([5]) which is 
distinct from the typical exponential decay apart from 
branch point singularities. Here it is important to note 
that this effect is not only observable exactly at the pa- 
rameters of the EP but can rather be seen in a large 



vicinity. In our study we found it still for a relative offset 
of the parameters of S = 10^^. 

We were furthermore able to demonstrate that it is 
possible to excite an adequate superposition of the eigen- 
vector at an exceptional point and its associate counter- 
part in a realistic physical situation such that the unique 
time signal becomes observable in atomic spectra. The 
quadratic term significantly influences the survival prob- 
ability we found for the hydrogen atom in crossed electric 
and magnetic fields. It is an effect that leaves clear sig- 
natures in the decay of resonances in quantum systems 
and obviously opens a new possibility to detect an un- 
ambiguous fingerprint of an exceptional point accessible 
with experimental methods. It might such facilitate the 
first experimental detection of exceptional points in a 
true quantum system. 

In this article the time evolution of the resonances ex- 
cited with the laser pulse (|14l) is evaluated with the sur- 
vival probability calculated from the spectrum via the 
c- product [2l|, [2I]. Both its modulus and phase have 
experimental consequences ^2^. In the realistic physi- 
cal situation shown in Fig. [3] the survival probability de- 
scribes the decay of a resonance state in time. Thus one 
has to measure the decaying occupation of such a state. 
This can presumably be detected with a second laser as 
already indicated 'l6'|. Furthermore, an extraction of the 
time signal from the spectrum is possible and has al- 
ready been used for microwave cavities (l9j via a Fourier 
transform. In this context we should mention that the 
ultra-strong magnetic and electric fields we have used in 
our calculations are due to our computational limitations 
to exceptional points associated with low-lying resonance 
positions. However, the phenomenon discussed here ap- 
pears for highly excited resonances as well. In such cases 
much weaker and feasible external fields are required to 
observe the unique survival probability of a wave packet, 
which initially populates mainly the two almost degener- 
ate states associated with the exceptional point. 
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